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Introduction 


The analog computer has proved to be a very valuable tool in 
solving engineering problems. Many of these computational problems 
require the use of trial-and-error or searching procedures if solutions 
are to be obtained. The use of conventional trial-and-error techniques 
by the computer operator is often quite time consuming, particularly 
if a large number of solutions is required. 

In recent years "hybrid" computers have been developed which 
employ both analog and digital elements. The searching operation re¬ 
quired in solving problems of the type mentioned above can often be 
performed automatically and at high speed on a hybrid computer. If 
the searching operation is to be implemented entirely by the computer, 
it is necessary to specify a sequence of operations which the computer 

A // 

is programmed to perform repeatedly. The terms iterative procedure 
and algorithm are used to denote this sequence of operations. 

Two-point boundary value problems belong to the class of prob¬ 
lems for which searching procedures are required. In this paper a 
simple example of a two-point boundary value problem is considered. 

An iterative procedure is devised for this example and a computer 
circuit for implementing this procedure is described. Numerical re¬ 
sults obtained by using this circuit are presented. This example pro¬ 
vides a convenient means for illustrating the mathematical difficulties 
which arise in the use of iterative procedures. 

The Example Problem 

For the one-dimensional flow of heat in a slab, the temperature 
distribution is determined by solving the equation. 


C(x) 
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+ f(x, t). 
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subject to appropriate initial and boundary conditions, where: 


0(x, t) = temperature, 

C(x) = specific heat, 

K(x) = thermal conductivity, 

x = distance along the slab, 0 s x s 1, 
t = time, and 

f(x, t) = rate of heat flow per unit volume into or out of 
the medium. 

Equation (1) is a linear partial differential equation. For the 
purposes of this paper it will be assumed that f(x, T ) s 0. The method 
of separation of variables will now be used to separate equation (1) 
into two ordinary, linear differential equations. In this method a solu¬ 
tion of the form 

0(x, t) = X(x)T(t) (2) 

is assumed. Substituting (2) into (1) and dividing the resulting equation 
by X(x)T(t) gives 

[T(t)] _ 1 -^3 = [C(x)X(x)] _1 ^ [K(x) ] . ( 3) 

The only way in which this equation can be satisfied as x and r vary is 
for each side of this equation to be equal to a constant, say a. It can 
easily be shown from physical considerations that a must be negative. 
Setting \ = -a (so that X > 0), the following equations are obtained 
from (3): 

5^+\t = o, ( 4) 

s[ K(x| f J + " c(x)x = 0 - < 5 > 

The solution of equation (4) is straightforward. However, 
equation (5) must be solved subject to the boundary conditions at 
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x = 0 and x = 1. For simplicity, assume that the boundary conditions 


are 

X(0) = X(l) = 0. 


( 6 ) 


The solutions of equation (5) will only satisfy the prescribed boundary 
conditions (6) for certain discrete values of say X , where n = 1, 2, • • •. 
The numbers X are called eigenvalues. Let X (x) denote a solution 
of equation (5) which is obtained using X = X^. Then X^fx) satisfies 
the boundary conditions (6). The function X n (x) is termed an eigen¬ 
function or normal mode. 


These concepts will now be illustrated by considering the case 
where the slab is homogeneous, i.e., where C(x) and K(x) are con¬ 
stant for OSxSl, Suppose that C(x) = K(x) = 1. For this case equa¬ 
tion (5) becomes 


and the solution of 


d 2 X 


+ XX = 0 


dx 

this equation is 


(7) 


X(x) = A sin '^/x'x + B cos N ^x'x. 

Since X(0) = 0, it follows that B = 0. The other boundary condition 

requires that A sin '^x' =0. In order to obtain a non-trivial solution, 

it is necessary that sin x^X = 0, or = nir. Hence the eigenvalues 

for this case are X = (mr) 2 , n = 1, 2, • • • . The function X (x) = 

n m 

A sin ^x m x is the mth normal mode with amplitude A. 

The problem to be solved is the following: Find the first N 
eigenvalues and corresponding normal modes of equation (5) subject 
to the boundary conditions given by (6). 

The Proposed Method of Solution 

Once the problem has been stated in mathematical terms, the 
next step is the formulation of a logical procedure to be programmed 
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on the computer. A procedure will now be devised for the problem 
described above. 


Consider the solution of equation (5) on the hybrid computer. 
The independent variable in the computer solution is computer time. 
The boundary condition X(0) =0 becomes an initial condition for the 
computer solution. The other boundary condition, X(l) = 0, must be 


satisfied by trial-and-error adjustment of the parameter X. Note 

dX 

, which is needed for the computer solution, 
x=a_ _ dX 


that the choice of 


dx 


is arbitrary, since the value of 
not the shape of the solution X(x). 


x=0 


affects only the amplitude and 


Suppose that the computer is seeking the eigenvalue X . The 

computer must use some criterion in correcting the parameter X on 

each iteration to insure that when |X(1) | < € (c small determines the 

accuracy to which the boundary condition is satisfied), X(x) is close 

to X (x) and not some other normal mode. The criterion to be used 
n 

is based on the zero-crossing property of the normal modes: namely, 
the larger the eigenvalue X^, the closer together lie the zeros of the 
normal mode, X^(x). This property will now be shown. * 

Let X . and X, be two distinct eigenvalues of equation (5) for 
1 

the boundary conditions (6); let X^(x) be a solution of equation (5) 

with X = Xand let X. (x) be a solution of equation (5) with X = X . 

j k k 

(Note: This proof is not restricted to the boundary conditions given 
by (6), but holds for most of the usual boundary conditions. ) Multi¬ 
plying the equation for X.(x) by X, (x) and that of X, (x) by X.(x) and 

] x x 3 

subtracting the second of these equations from the first yields 
dX. dX 

V K(X > + b CW Vk - X j [K(x) “dx 1 - \ C(x)X j X k = °- 


*For a more detailed discussion of this and related properties, see 
reference [1]. 
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or 


, dX. dX 

d - [KM (X k - X. —)] = (X k - X.) C(x) X.X k . (8) 

For some range of values of X the solution, X(x), of equation (5) goes 

through zero at least once for 0 < x < 1. Let X lie in this range, and 

let a be the smallest value of x > 0 such that X.(x) = 0. Clearly a < 1. 

J 

Integrating equation (8) from x = 0 to x = a and applying the conditions 

X.(0) = X, (0) = X.(a) =0 gives the equation 
1 K J 


dX. 


x k KW ^ 


x=a 


= (\ -X.) \ C(o-)X.(<r)X. (cr)do- 
k J <V J k 


(9) 


dX. 

Suppose that X.(x) > 0 for 0 £ x £ a. Then < 0. Also suppose 

that X k (x) does not change sign for 0 £ x £ a. For the sake of argu¬ 
ment let X k (x) be positive for 0 £ x £ a. Now K(x) > 0 and C(x) > 0 for 
0 £ x £ 1. Inspection of equation (9) shows that, under the above hypo¬ 
thesis, the right hand side is positive while the left hand side is nega¬ 
tive. This is impossible. Hence X fc (x) must change sign at least once 
for 0 £ x £ a. This same argument can be used to show that X (x) 

JK 

must change sign between any two successive zero-crossings of X.(x). 


Let the eigenvalues be ordered numerically so that 0 < Xj < X 2 

<\ 3 <. . • . (This ordering is possible, since the eigenvalues are all real.) 

Then, if X (x) and X (x) are solutions of equation (5) with X = X 
n n+i n 

and X = k n+1 respectively, it follows from the above argument that 
X n+j (x) has at least one more zero-crossing than X^(x) for 0 £ x £ 1. 
Since the solution, X(x), of equation (5) depends continuously on X, it 
follows that X n+1 (x) can have at most one more zero-crossing than 
X n ( x ). If this were not the case, then there would be a value of X, 
say X., such that X^ < X < X^ +j and the solution of equation (5) with 
X = X would satisfy the boundary conditions (6), i. e., X would be an 
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are succes- 


eigenvalue. But this contradicts the fact that X and X 

n n+i 

sive eivenvalues. Hence X n+i (x) has just one more zero-crossing 
than X n <x) for Oi xSl. It follows from the definition of an eigenvalue 
and the fact that the eigenvalues are ordered, that X n (x) has n zero- 
crossings for 0 < x £ 1. (Note that X (0) = 0 is not considered as a 


n 


zero-crossing of X (x).) 


n 


Suppose that the nth eigenvalue, X , and a corresponding nor 


mal mode, X (x), are to be determined. Let X 1 denote the ith approxi- 
n n 

mation to X ; let X (x) be a solution of equation (5) obtained with 


X = X X n ; and let x^ denote the value of x for which the nth zero-crossing 

of X^(x) occurs (i. e., let x^ be the nth positive root of the equation 

X 1 (x) = 0). Note that if X* = X , then x 1 = 1. It follows from the dis- 
n n n n 

cussion of the zero-crossing property of the normal modes that, if 

x^ > 1, then X^ < X^ and hence X^ (i. e., the (i + 1) approximation to 

X ) should be chosen such that X i+1 > X* . Likewise, if x 1 < 1, then 
n i+1 n n . n 

X* > X and X should be chosen so that X* +1 < X 1 . This suggests 
n n n ^ n n &e> 

that one way of choosing X^ is by means of the recurrence relation 


. 1+1 = X 1 + k(x 1 - 1) , 
n n n 


( 10 ) 


where k > 0 is a constant. 


The problem proposed earlier is that of determining the first 
N eigenvalues and corresponding normal modes. The iterative pro¬ 
cedure to be considered is based on equation (10). In this procedure 
the computer must perform the following sequence of operations: 

1) Store the number n in a reference counter; 

2) If n > N, stop the program; otherwise proceed to step 3); 


3) Solve equation (5) with X = X 1 for X^(x) subject to the 

n ^ivi 

= A where A / 0; 


dX 1 

initial conditions X„(0) = 0 and —— 

n dx 


x=0 


4) Count the zero-crossings of X 1 (x) and compare this number 


n 
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to the number n stored in the reference counter to deter- 
i 

mine x ; 
n 

5) If |x - 11 > e (where « > 0 is specified a priori), use 
equation (10) to obtain X^ + *and repeat steps 3), 4), and 5.); 

6) If |x* - 11 £ €, record X 1 and X 1 (x) for 0 £ x £ 1, add one 

n n n 

to the number stored in the reference counter, and return 
to step 2). 

Convergence of the Iterative Procedure 


The iterative procedure for determining X n generates a se¬ 
quence (X^, X*, . . ., X n ) of approximations to X . An important 
question to be considered is whether or not this sequence actually 
converges to X . Useful information can be gained by linearizing 

equation (10) about X . 

n 

Let X(x) be a solution of equation (5) for some positive value 

of the parameter X. Suppose that the equation X( x) =0 has at least 

m positive roots for 0 < x < oo where m > 1. Let x. denote the jth 

positive root of X(x) = 0 where 1 £ j £ m. It was pointed out earlier 

that X(x) depends continuously on X. It follows from this fact that 

x = x.(X) is a continuous function of X. This is true for each i. 

J J i i J 

Hence x = x (X ) and x (X ) = 1. 
n n' n' n v n' 

Suppose that x^(X) has a Taylor's series expansion about the 

point X = X . Then 
n 


x (X) = 1 + 
n 


dx 

n 

"dX~ 


(X - X ) + r(X - X ), 
n n 

n 


(ID 


where (X - X^) r(X - X^) -*• 0 as \ -*• X . Substituting equation (11) 
with X = X^ into equation (10) gives 


i+i 


dx 


n 


X +k 
n 


n 


dX 


X -X 


(x 1 -X ) +kr(x x -X ) , 
n n n n 


( 12 ) 


n 
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so that the linearized difference equation corresponding to equation 
(10) is 


6X 


i+l 

n 


dx 


1 + k 


n 


dx 


x=x 


n. 


6X 


n 


where 6X. 1 = X 1 - X . The solution of equation (13) is 
n n n 


6X 


n 


dx 

1 + k dV 1 


II 


n _ 


6X' 


n 


(13) 


(14) 


where X” is the initial approximation to X^. 

It can be shown that, if 6X° is sufficiently small, then the 
solution of equation (10) (which has X° as its initial value) converges 
to zero provided that all solutions of equation (13) converge to zero. 
It is clear from equation (14) that solutions of equation (13) converge 



dx 


dx 

to zero if and only if 

1 + k dx" 

X =x 

n 

< 1. Since 

dX 


< 0 , 


X =X 


n 


this condition requires that 
0 < k < 2 


dx 

n 


dX 

X =x 

n 


(15) 


Thus, if the value of k chosen is too large, the iterative procedure 

will not converge. On the other hand, if k is chosen conservatively 

small, the convergence may be quite slow. Furthermore, since 
dx 


n 

dx" 


X =x 


depends on n, the range of k for which solutions of equa- 


n 

tion (13) converge may vary significantly for different n. 


To clarify these points, the case of the homogeneous slab 
described by equation (7) will be considered briefly. It follows from 
the earlier di scussion of this case that x is determined from the 


equation 


n 
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x = nir. 
n 


_! 

Therefore, x^(\) = mr(\) 2 and 


dx 

n 

d\ 


\=X 


= -imr(\ n ) 


1 

2 = 


n 

-j(n-rr) , since = (nir) 2 . For this case, condition (15) takes the form 
0 < k < 4(mr) 2 . Suppose that k = 2it 2 . Then for n = 1, equation (14) 
shows that 6\ l = 0 for any 6X.j. On the other hand, suppose k = 2tt 2 and 
n > 1. Then equation (14) becomes 


5X 


i 

n 


(i - 4) l6x ° 
n 2 n 


and it is evident that the convergence is quite slow for large n. 

This problem of selecting a useful value of k is not quite as 
serious from a practical standpoint as the above discussion would 
seem to indicate. The reason is that the computer can usually be 
programmed to automatically increase or decrease the value of k (at 
least within some range of values). Alternatively, the computer oper¬ 
ator can make appropriate changes in k whenever such changes are re¬ 
quired. 


Another question which presents itself is whether or not the 
boundary condition at x = 1 is satisfied with sufficient accuracy when 
l<-l|s £ for a specified «. It may well be that X^(l) is relatively 
large even though |x n - l| is quite small. A better iterative pro¬ 
cedure for this problem would be one based on reducing |x^(l) | di¬ 
rectly. However, the computer circuit required for this type of 
algorithm is more complex than the circuit required for the iterative 
procedure described earlier. 


The Hybrid Co mputer 

Now that an iterative procedure has been formulated for the 
problem, the next step is to obtain a computer circuit for implement¬ 
ing this procedure. First, however, the hybrid computer will be dis¬ 
cussed briefly and the symbols required for programming the iterative 
procedure will be described. 
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The elements in a hybrid computer can be categorized as analog, digi¬ 
tal, and analog-digital, depending on the class of signals with which 
the elements deal. For the problem considered in this paper only the 
patchable digital elements to be described below are required. How¬ 
ever, in problems which require numerical calculations and complex 
schemes for control of the analog-digital elements, a stored program 
digital computer with flexible input-output routines is needed. 

The symbols used for programming the iterative procedure 
are shown in figure 1. Light-weight lines and lower-case letters are 
used to designate the analog signals, while heavy-weight lines and 
upper-case letters are used for logic signals. A logic signal. A, is 
either in the state A = 1 or in the state A = 0. A is used to denote the 
complement of A, i. e., if A = 1, then A = 0 and if A = 0, A - 1. 

The operation of the elements depicted in figure 1 will now be 
described. 

a) Summer: +eo = <*(10ei + ez + i). The standard feedback 
is 0.1MS2 (i. e. , a = 0.1). Note that both plus and minus 
output signals are available. This type of element is 
called bi-polar. 

b) Potentiometer ! eo = kei. 

, -i 

c) Multiplier: l = 10 xy. 

d) Integrator : Each integrator is separately controlled by the 
logic signals A, B, C, and D. Electronic switches are used 
for high speed mode changing. 

Reset Mode: A = 0, B = 0, 1; eo - e. 



o 


Hold Mode: A = 1, B = 1. 

Time Scale! C = 0, D = 0; real time (integrator gain = a). 


C = 1, D = 0; speed up computer time by factor 
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of 10 

C = 0, D = 1; speed up computer time by factor 
of 100. 

The standard feedback is a 1 (ofd. capacitor 
(i. e., a = 1). 

e) Switch (Single pole, double throw): 

A = 0; i = 10ei 
A = 1; i = 10e 2 

f) Comparator ! 

B = 0} A = 1 if ej + e 2 > €; 

A = 0 if ej + e 2 < -€; 

A = 0,1 if | e 2 + e 2 | < €, i. e. , for | e 2 + e 2 | < €, 
the output A depends on the past history of 
ex + e 2 so that A has the value it last had when 
l e i + e 2 | > c. (The small dead space e and 
associated hysteresis loop provides noise sup¬ 
pression). 

B = 1: Locks A at value present when B = 0 -* 1 where 
0 -*• 1 denotes a transition from logic 0 to logic 1. 

g) Track-Transfer (Memory element): 

A = 0, B = 0, 1: eo = e (Initial condition mode) 

A = 1, B = 0, 1: hold eo constant 

A = 1, B = 0 —1: transfer (ej + e 2 ) to e 0 and hold e 0 
constant until next B = 0 -*1. 

h) Input Button : 

B = 0: A = 1 if button down, 

A = 0 if button up. 

B = 1: A = 0, 1, push button to reverse A. 
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i) Pulser (Variable period): 

A = 1, 0; B = 0 

A = 0 -* 1: B = 1 for T seconds 
C = 0, D = 0: T = 1 second 
C = 1, D = 0; T=10 milliseconds 

C = 0, D = 1: T = 100 milliseconds 

j) Pulser (Fixed period): 

k) Flip-flop : 

A = 0, B = 0: D = 0,1; 

A = 1, B = 0: D = 1 
A =0, B = 1: D = 0 
A = 1, B = 1 not used 
A = 0, B = 0, C = 0-*1: reverse D. 

l) Or Gate : 

D = 0 only if A = 0, 6 = 0, and C = 0. 

The Computer Program 

A block diagram of the proposed computer program is shown in 
figure 2. The circuits for the various blocks are contained in figures 
3-8. These circuits are described below. Timing diagrams are in¬ 
cluded with most of these circuits. These timing diagrams are not only 
helpful in explaining the behavior of a circuit, but are also useful in de¬ 
signing a computer circuit to perform a specified operation. 

In the computer solution of this problem, the independent variable 
is computer time t. The relationship between t and x is t = tx so that 
t =~t corresponds to x =1. Equation (5) can be written in terms of t as 



( 16 ) 


where y(t) = X(t 1 t). 

The on-off control for the program is pushbutton 1. When this 
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control is in the "off" state (i.e., S = 1) the integrators are in reset, 
the track-transfers are in the initial condition mode and the flip-flops 
are cleared. 

The timing circuit shown in figure 3 is the clock which controls 
the reset and operate periods. It also furnishes the independent variable 
t for use elsewhere in the circuit. When B = 1, the integrator goes into 
the operate mode and generates the ramp ct. For ct < 100 volts, the 
outputs of Cl and PI are each logic zero. When ct exceeds 100 volts, 
the output of Cl changes from 0 to 1 and triggers PI. This in turn puts 
the integrator in the reset mode. This clock can be interrupted by- 
placing the integrator in the hold mode. This is done by signal D when¬ 
ever a zero-crossing of y*(t) occurs (where y^(t) is a solution of equa¬ 
tion (16) obtained for X = X^). The signal F comes from the manual and 
programmed stop circuit and F = 0 unless a stop has been requested. 

The relationship between the various signals is shown in the 
timing diagram of figure 3. 

Figure 4 sho-tfrs the zero-crossing counter. The reference 
counter is TT2. The output of TT2 is 10(n - jr) volts where n is the 
number of zero-crossings of y (t) for 0 < t < t. The initial condition 
of -5 volts on TT2 is used to insure an accurate switching of C3 when 
the nth zero-crossing of y^(t) occurs. Once X^ and y (t) have been de¬ 
termined and recorded, the print out routine causes G = 1 for a short 
period of time. The signal G = 0 - 1 triggers TT2 and its output changes 
to 10(n + j) volts. 

The number of zero-crossings of y^(t) is counted by TT1. On 
each zero-crossing, the output of C2 changes state and triggers either 
P2 or P3. Thus for each zero-crossing, D = 1 for a period of several 
milliseconds. The signal D = 0 -* 1 causes 10 volts to be added to the 
output of TT1. The signal D = 1 is also used to put the clock (figure 3) 
and the integrators used in the simulation of equation (16) in the hold 


- 13 - 




mode while TT1 is being updated. This insures an accurate determin¬ 
ation of (t* - t), if ct 1 < 100 volts (where t 1 = t x 1 ). The output of TT1 
n n n n 

is reset to zero during each reset period. 

At time t = t 1 (if ct 1 < 100 volts) the output of C3 changes state, 
n n 

i.e., A = 0 — 1. This signal is used to initiate the parameter correc¬ 
tion routine. 

The parameter correction circuit is shown in figure 5. When 
|t^ - t| < e, the signal H is logic one and H is used to prevent any 
change in x\ However, if |t n " t | - « as is generally the case, H = 0 
and the value of X is updated when A s 0->lorC = 0-*l. If ct^ <100 volts, 
A = 0 -*• 1 when t = t^ and this signal causes X to be updated according to 
the equation 

X i+1 = X* + kc(t* - t) . (17) 

n n n 

On the other hand, if ct 1 2: 100 volts, A = 0 for 0 £ ct s 100 volts and 

n 

C = 0 -* 1 causes X to be updated as follows: 

X 1+1 = X 1 + k(100 - ct) . (18) 

n n 

The signal E comes from the automatic scale changing circuit. 
When E = 1 the gain constant k in equations \17) and (18) is replaced by 
0.1k. 

Since the problem stated that the first N normal modes are to 
be determined, the computer must be programmed to stop after Nth 
mode has been recorded. It is also convenient to have a manual stop 
button. This, when depressed, causes the computer to stop after it 
has determined and recorded the normal mode that it is currently seek¬ 
ing. This manual stop is useful during check-out for example. The 
circuit for performing the stop operation is shown in figure 6. 

The end point error test circuit is shown in figure 7. The 
comparator C5 compares £ |t - 11 to e . If ct^ < 100 volts, A = 0 — 1 
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at t = t 1 and the signal A = 1 is used to latch C5. If ^It 1 - t| < e, the 

n n 

flip-flop is set, i. e. , H = 0 -*»1, and the print out routine is initiated. 
For simplicity the details of this operation have been omitted. At the 
termination of the print out procedure G - 1 for a short period. This 
signal is used to set H - 0 and update the reference counter, TT2. 

The signal C = 1 is used to latch C5 during the reset period to insure 
that false triggering of the flip-flop does not occur. 

Figure 8 shows the automatic scale changing circuit. The 
addition of this circuit allows a significant increase in the number of 
eigenvalues that can be determined before an overload occurs. 

In figure 8 the integrator is used as a track-hold device. It 
is normally in the track mode. When the output of TT3 exceeds 100 
volts, J = 0 -* 1. This triggers the pulser which latches C6. The sig¬ 
nal J = 1 also sets the flip-flop, puts the integrator in the hold mode, 
and puts TT3 in the initial condition mode. The output of TT3 is now 
set to one-tenth of the value it formerly held. The signal E = 1 is 
used to decrease by a factor of ten the size of the integrating capacitors 
in the circuit which solves equation (16) and to decrease the gain in the 
parameter correction circuit by a factor of ten. The signals A and C 
are used to latch C6 so that the scale changing does not occur while 
TT3 is in the transfer mode. Additional circuitry is required if it 
is desired to make another scale change when and if the output of TT3 
exceeds 100 volts for the second time. 

Numerical Results 


Numerical results were obtained for the case of the homogene¬ 


ous slab described by equation (7). For this problem equation (16) re¬ 
duces to 



+ t 2 \y = 0. 


(19) 


The circuit for solving this equation is shown in figure 9. To obtain 
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isu 
N _ . X 


used as the parameter input 


good amplitude scaling, 100 t 

signal to each of the multipliers. Let y 1 = t 

n M n 

lem it is convenient to modify the iterative procedure slightly. The 



For this prob- 


computer was programmed to seek y^ rather than t 
tive procedure used was based on the equation 


n 


The itera- 


i+i _ i , , ,.i 7, 

y = y^ + kc(t^ - t). 


n 


n 


n 


( 20 ) 


It is helpful at this point to linearize equation (20) about y . 
Since the derivation of this linearized equation parallels the derivation 
of equation (13), the details are omitted. The desired linearized equa¬ 
tion is 


Sy 


i+i 

n 


dt 


1 + kc —^ 


ay 



Y = v 

1 ’n 


6y 


n 


( 21 ) 


where 6y 


n 


y - y . 
*n T n 


0 < k < 2 


dt 


n 


C dy 


y-y 

1 l n 


All solutions of this equation converge to zero if 

_1 -i dt n 

. Since t n (y) = mr(y) , = -n-rrfyj 


y=y 


n 


= -t 2 (nir) 1 (using the fact that y n = rnr(t) *). Thus the above condition 
on k can be written as 0 < k < 2 | -ct 2 (nir) 1 1 1 . In obtaining the numeri¬ 
cal results the values c = 6. 667 volts/sec and t = 10 sec. were used. For 
these values of c and t and for n = 1, it is necessary that 0 <k < 2mr(ct 2 ) 1 
= 0.00942 if all the solutions of equation (21) are to converge to zero. 

The value of k actually used was k = 0.006. 


dt 

Substituting the above values for c, k, and —^ 


tion (21 ) yields 


Y = Y 


n 


6y^ +1 = [1 - 1.274(n) ‘*]6y* , 


into equa- 


( 22 ) 


and the solution of this equation for the initial value 6y^ is 

-l _i 


Sy^ = [1 - 1. 274(n) ~ f 6y° . 
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This shows that the number of iterations required to obtain a solution 
of specified accuracy will increase as n increases. 

The numerical results obtained for this problem are presented 
in Table 1. Inspection of the data shows that equation (22) gives a 
good qualitative indication of the type of behavior to be expected for 
in the vicinity of y . 

Conclusions 

In engineering applications computational problems often arise 
which require the use of trial-and-error procedures if solutions are 
to be obtained. For many problems of this type conventional trial- 
and-error methods can be replaced by iterative routines which can 
be implemented on a hybrid computer system. 

In this paper a simple example is considered. An iterative 
procedure is formulated and the mathematical difficulties which arise 
in the use of this procedure are discussed. A computer circuit for 
implementing this iterative routine is described and some numerical 
results obtained by using this circuit are presented. 

Reference 1 

P.M. Morse and H. Feshback, "Methods of Theoretical Physics, 

Parti, pp 719-7 24, McGraw-Hill Co., Inc., 1953. 
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Figure 1 
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Figure 2. Block Diagram 
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Figure 3. Timing Circuit 
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Figure 5. Parameter Correction Circuit 
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Figure 6. Manual and Programmed Stop Circuit 
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Figure 7. End Point Error Test Circuit 
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Figure 9. Uniform Slab Simulation 














































Table 1. Numerical Results for Homogeneous Slab 


V! = .3142 

y 2 = .6283 

Ys = -9425 

-1.2566 

100 vj 

100 y\ 

100 Y 3 

100 Y4 

10.. 00 

31.41 

62.81 

94. 24 

29.97 

51.38 

82.79 

11.42* 

31.89 

60. 29 

88.32 

11.82 

31.29 

61. 97 

90.99 

12. 07 

31.44 

62. 52 

92.42 

12. 24 

3 1.40 

62. 71 

93.21 

12. 34 

31.41 

62.78 

93.66 

12. 41 


62. 80 

93.91 

12. 46 


62. 81 

94.05 

12. 49 



94. 13 

12. 52 



94. 18 

12. 53 



94.21 

12. 54 



94. 23 

12. 55 



94. 24 

i 



*Rescaled to 10 . 
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